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Abstract 

The series solution of the radial part of the Schrodinger equation for simultaneous coulomb and 
harmonic potential involves three-term recursion relation and is thus difficult to solve for bound 
states. We have suggested a simple method to solve for low lying states. Finite polynomial solutions 
exist only if the coulomb and oscillator potentials are nontrivially related. 
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I. INTRODUCTION 

In the quantum mechanics text books [1], the radial part of the Schrodinger equation for 
a central potential problem is solved by Frobenius' method where in general a the two-term 
recursion relation relates one coefficient of the series to another one. Normally one needs to 
truncate the series to a finite polynomial to have normalizable bound state wavef unctions. We 
also get the energy eigenvalues from the condition of the series truncation. But we find very rare 
comments about the recursion relation involving more than two coefficients as they are difficult 
to solve. In the a popular quantum mechanics book, in the context of hydrogen atom problem, 
the three-term recursion relation has been commented as "enormously more difficult to work 
with" compared to a two-term recursion relation[2]. We can avoid the three-term recursion 
relation in hydrogen atom by choosing the functional form of the radial wavefunction from its 
behavior at small and large r. Here we discuss one such interesting example where we cannot 
avoid a three-term recursion relation. Consider the Schrodinger equation in three dimensions 
when both coulomb and harmonic oscillator potentials are present together, i.e., the potential 
is given by 

V(r) = —a — h -mu 2 r 2 (1) 

(where a = e 2 /(47reo) for hydrogen atom problem). The series solution of the radial part 
for this potential involves a three-term recursion relation. For the bound state solution the 
series needs to be truncated to a polynomial. The main difficulty of this problem is that the 
analytic solution for any arbitrary energy level in general form is very difficult to obtain. In j^] , 
a similar type of problem was addressed with anisotropic harmonic oscillator with frequency 
ratio of the oscillators in different directions chosen in such a way that separation of variables 
in the parabolic co-ordinates works. If one takes isotropic oscillator, then it is not possible to 
make the separation of variables for all three parabolic coordinates as done in The authors 
of that paper tried a series solution which again resulted in the three-term recursion relation 



which cannot be solved analytically and they solved numerical 



case when the coulomb term is zero. Hall, Saad and Sen {4] solved the isotropic case in d- 
dimensions with more rigorous mathematical approach. They showed that to have polynomial 
solutions the parameters in the potential need to satisfy specific conditions depending on the 



y in two dimensions for a special 



order of the polynomial. 

In this paper, we have suggested a very simplistic method to obtain the low lying energy 
eigenvalues without much mathematical complications. Converting the recursion relations in- 
volving three coefficients to relations involving only two, we can find out the conditions for 
the series to be terminated to give normalizable bound state wavefunctions and the energy 
eigenvalues. 



II. THE RADIAL EQUATION 

We write the Schrdinger equation in spherical polar coordinates. After the separation of 
variables, the angular part of the Schrodinger equation can very easily be solved and can 
be found in any standard quantum mechanics book 1|, the solution is given by the spherical 
harmonics Yi m {6, 0). The radial equation is given by 

h 2 d 2 u(r) r a 1 , 9 h 2 1(1 + 1)-, . . . > . . 

where I is the azimuthal quantum number. The effective potential 



a 1 o o h 2 1(1 + 1) 



Let us define 
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*=-jj-> P = kr, Po = W ^d Pl = — . (4) 



Then the radial equation (Eqj2]) can be written as 

d 2 U [22 1(1 + 1) Po 



dp 2 ^ p+ — 2 — -p- ^ (5) 

The asymptotic behavior is determined by 

CiZ)_i/2(\/2pi p) + C , 2-D-i/2(«v / 2pi P ) (where Ci and C2 are constants). With the condition 
that the wavefunction vanishes at infinity, the solution for p — > 00 goes as 

-PiP 2 /2 



u ~ e~ Plp 11 (7) 



whereas the bevavior at small distances (p — > 0) is given by 

u ~ p l+1 (8) 
So, we assume that the general form of radial wavefunction is given by 

u(p) = p l+1 e- plp2/2 v(p) (9) 

With this substitution, Eq.(j5]) becomes 

d^i) civ 
p— + 2(1 + 1- Pip 2 )j- p + [po - p{-l + (21 + 3) Pl }]v = (10) 

A. series solution 

Substituting 

oo 

v( P ) = j2 c *p t (ii) 

in Eq. ffTOl) we get 

oo oo oo 

+ l)c m p* + 2(Z + 1) + 1)^+!^ - 2 Pl ^ z Q p m 

oo oo 

+ po Yl C ^ ~ t( 2/ + 3 )-°i - !] £ C ^ i+1 = ( 12 ) 

i=0 i=0 

The coeffients q satisifies the recursion relation 

[q + 2pi(z- l)]cj_i - ppCj 

{% + !){% + 2 + i) {V6) 



where a = — 1 + (21 + 3)pi and 



Ci = ~2(rh) Co - (14) 



Looking at the solution for large i, the approximate beviour of the series solution is given by 

OO rfz OO 1 



Each of the terms in the above expression goes as e pip and thus spoils the asymptotic behavior 
of the radial wavefunction u(r). So, we need to truncate the series. Note that here the recursion 
relation Eq.( fl3l) involves three coefficients Cj+i, q, q_i and the truncation is very tricky. If 
we want to truncate the series in the conventional way that is if we set the coefficient c n+ i = 
it does not guarantee the termination of the series. On top of that, if we set c n+ i = then it 
implies from Eq. (fl3|) that 

c n o + 2/9i (n-1) , . 
= (16) 

C n -1 PO 

increases as n increases i.e., the series does not converge which contradicts the original recursion 
relation! In its present form, we cannot set Eq. (1131) to be zero. 

To truncate the series we need the recursion relation with only two coefficients. For this 
purpose we rewrite the first few coefficients in term of the lowest order coefficient cq\ 

Po 

ac - pqCi (21 + 2) a + pi 
° 2 ~ 2(2/ + 3) ~ 2(2/ + 2)(2/ + 3)° 

(a + 2pi)ci - ppc 2 p [2(2/ + 3)(a + 2pQ + (2(1 + I) a + pi)} 
° 3 3(2/ + 4) 6(2/ + 2)(2/ + 3)(2/ + 4) C ° [ ' 

etc. Alternatively, one can also write down the coefficient c n in terms of c n _i as 

Po 

Cl " ~2(7Tl) Co 

ac - p ci 2(1 + l)a + pi 
° 2 ~ 2(2/ + 3) ~ 2(2/ + 3)p Cl 

(a + 2 Pl ) Cl - p c 2 _ p [2(2/ + 3)(a + 2 Pl ) + (2(1 + l)a + gg)j 
C3 3(2/ + 4) (2(/ + l)a + p2) 3 (2/ + 4) ° 2 1 j 

and so on. It is not possible to write down the recursion relation for a general term c n in these 

forms, but we can write down as many terms as we wish. So, in place of a general solution 

let us look at the low lying solutions. The lowest possible term that can be set to zero is c 2 

(this corresponds to % = 1 in Eq. (ll3p ). If c 2 — then C3 = only if the coefficient of C\ in the 

recursion relation for C3 i.e., (a + 2p\) = which puts an additional constraint on the energy 

eigenvalues. The conditions that c 2 = is given by 

2(1 + l)a + pl = 0. (19) 



Since Pi — ^ and pi = t§§t, Eq. (jT9|) gives the energy eigenvalues 

Ea = ( 2l + ,)-^ + ww . (20) 
The condition that C3 also be zero is given by 

a + 2 Pl = (21) 
which gives another expression for energy eigenvalue 

£i,= ~(2Z + 5)/iu;. (22) 

Eq. pTI) together with Eq.(IT9l) gives the condition 

2 

/? = = / + 1 = integer (23) 

Thus, the infinite series can be terminated into a finite polynomial only if the parameters in 
the effective potential satisfy a nontrivial condition! The above results exactly agree with the 



results in (4]. The minimum value of the ratio (3 is one thus this method is not applicable with 
a = 0. With this condition the energy eigenvalue Eq. (l20|) becomes 

^ /, 5.^ 3^ ma 2 , , 

E u = (l + -)?iLj = -HLJ + 1 p- (24) 

It is interesting to note that though we have never assumed any particular value for /, the final 
expression for the energy (Eq.(l2l])) corresponds to I = in Eq.( l20|) . Now, let us set c 3 = 
(i.e., % = 2 in Eq. (fT3|) ) with the condition that a + 4pi = so that C4 also becomes zero and 
the series terminates to a second order polynomial. Then from Eqs.( fT7|) we get 



2{a + 2pi)(2Z + 3) + [2(1 + l)a + p 2 ] = (25) 

which gives 

_ 3(2/ + 3)(/ + 2) ma 2 

E21 ~ 2(3/ + 4) ^ + (3/ + 4)h 2 (26) 

and the condition on the oscillator frequency and the strength of the coulomb potential is given 
by the dimensionless ratio 

P = 41 + 5 (27) 
6 
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in agreement with |4|. With this condition the energy eigenvalue reduces to 

^ /, 7 Xa _ 9 ma 2 . . 

E 2l = (l + -)huj = -hw + — (28) 

which again corresponds to I = in Eq. fl26l) . If we set the coefficient c n = then the series 
terminates if the pre-factor of c n _i in the recursion relation for c n+ i is zero i.e., [a+2pi(n — 1)] = 
which gives the energy eigenvaule 

E nl = {n + l + -)Hu, n = l,2,3---. (29) 

This condition ensures that the ratio c n /c n _i in Eq.( !T6|) is zero and the divergence of the series 
does not arise. As the general condition for c n = can not be written down, we cannot 
determine the condition amongst u, a and I for arbitrary n . 



III. CONCLUSION 



We have proposed a very simplistic method to extract the low lying eigenvalues by truncating 
a series involving three-term recursion relation. We have considered the quantum mechanical 
problem with both harmonic oscillator and coulomb potentials. The series solution of the radial 
equation results in a recursion relation involving three coefficients and is difficult to truncate the 
series into a polynomial. We have shown that the low lying eigenvalues and eigenfunctions can 
be obtained by a very simple method by writing the coefficient c n in terms of the lowest order 
coefficient Cq or the preceding coefficient c ra _i. The series truncation needs an extra condition 
on the harmonic oscillator frequency and the coulomb strength. Depending on the order of the 
polynomial, the relation between u and a comes out to be different, but the dimensionless ratio 
f3 always takes only integer values. The energy eigenvalues can be written as purely harmonic 
oscillator or hydrogenic energy levels and the eigenvalues then depends on angular momentum 
I, on the other hand the energies can also be written in combination of both harmonic oscillator 
and hydrogen atom energies by completely eliminating I which is exactly same as the eigenvalue 
for I = without the additional constraint between a and u. Using our method we can obtain 
a general expression for the energy eigenvalue for terminating the series at any arbitrary n-th 
order polynomial, but we cannot evaluate the relation between a and u in the general form. 
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Thus this method provides a very simple way to extract the low lying eigenvalues when the 
recursion relation involves more than two coefficients. 
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